Quantum Treatment of Continuum Electrons in the Fields of Moving Charges 
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An ab initio, three-dimensional quantum mechanical calculation has been performed for the time- 
evolution of continuum electrons in the fields of moving charges. Here the essential singularity 
associated with the diverging phase factor in the continuum wave function is identified and removed 
analytically. As a result, the continuum components of the regularized wave function are slowly 
varying with time. Therefore, one can propagate continuum electrons to asymptotically large times 
and obtain numerically stable, well-converged ejected electron momentum spectra with very low 
numerical noise. As a consequence, our approach resolves outstanding controversies concerning 
structures in electron momentum distributions. The main conclusions are general and are illus- 
trated here for ionization of atomic hydrogen by proton impact. Our results show that in order to 
obtain correct long-time free-particle propagation, the essential singularity identified here should be 
removed from the continuum components of solutions to the time-dependent Schrodinger equation. 

PACS numbers: 34.70.+e, 34.50.Pi 



Tracing the evolution of continuum electrons in the 
presence of time-dependent external fields as they move 
from microscopic to macroscopic distances has been of 
great fundamental interest since the beginning of quan- 
tum mechanics [1 . Observations are made at distances 
where they can be detected by devices such as photon 
or particle detectors which also have macroscopic di- 
mensions. Substantial progress has been made in the 
measurements of continuum electrons owing to the de- 
velopment of the cold target recoil ion momentum spec- 
troscopy technique [2J [3J [U [5] , ^ or exam pl e - We describe 
a corresponding theoretical advance in this Letter. This 
advance is made possible by combining two lines of re- 
search, namely, the lattice time-dependent Schrodinger 
equation (LTDSE) method of Schultz and coworkers 
[5J [7J [B] and the removal of essential singularities of con- 
tinuum wave packets described in this Letter but implicit 
in the dynamical Sturmian theory of Ovchinnikov and 
coworkers [9]. 

Although standard numerical techniques successfully 
apply for time-dependent Schrodinger equation (TDSE) 
of bound states, they meet with many difficulties for the 
continuum components due to highly oscillatory phase 
factors exp[ir 2 /2t] as is evident from the free-particle 
propagator (see Eq.(2.5.16) in Ref.[T]). Because these 
phase factors have essential singularities for large dis- 
tances, both fundamental and numerical difficulties arise. 
Even with today's most powerful computers, propagation 
to large times without accurate treatment of the essential 
singularity is not possible. Here we present a technique to 
remove the singularity and apply it to a prototypical sys- 
tem treated in an explicitly time-dependent approach. To 
the best of our knowledge, no standard numerical meth- 
ods exist that can rigorously compute the time-evolution 
of quantum systems from atomic to macroscopic scales. 



With removal of the essential singularity, accurate nu- 
merical solution of the TDSE to asymptotically large 
times becomes possible. This is important for physi- 
cal processes where electron motion at large times, in 
contrast to just near the parent nuclei, plays a decisive 
role. For example, in ion-atom collisions, classical tra- 
jectory Monte Carlo calculations [10] show that cusps in 
the ejected electron spectra do not emerge until the tar- 
get and projectile species are separated by more than 
5,000 a.u., many times larger than the initial atomic di- 
mensions. Moreover, cusps are known on fundamental 
grounds to occur and are prominent features of all ioniza- 
tion processes where nuclear charges are not completely 
screened in the final state. For ionization by strong laser 
fields, rescattering of electrons from atomic or molecular 
species is important for multiple ionization and high har- 
monic generation [TT]. Such rescattering occurs due to 
oscillations of the laser field so that the electron periodi- 
cally turns around at large distances and recollides with 
the core of the species from which it originated. 

Thus, the key innovation illustrated in the present 
work is of far-reaching importance because the motion of 
continuum electrons in external fields is central to many 
processes which have only been treated classically until 
now. Specifically, we show how to treat these processes 
on an ab initio, completely quantum mechanical basis. 
To begin with, the electronic wave function is written in 
the form (atomic units are used throughout the paper 
unless otherwise indicated) 

*(r,t) = x(g,0)$(9,0), (1) 

where q=(rv sin 0)/f2 and t=— (£1 cot 9) / v 2 , varies from 
to 7T, and v and £1 are arbitrary parameters in general 
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but later v will represent a physical velocity. The factor 
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expH(rVcot6>)/2] (2) 



is highly oscillatory and has an essential singularity at 
l/g=0 (i.e. q — > oo). 

The regularized wave function ^(q, 0) however, has 
slow-varying continuum components and satisfies a reg- 
ularized time-dependent Schrodinger equation (RTDSE) 
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= «n*^, (3) 



where V(q,8) is an external, possibly time-dependent, 
potential. Time-dependent scale-transformations have 
been used to solve a model problem exactly [12] and a 
one dimension problem numerically |13j . but the method 
introduced here enables efficient and full numerical cal- 
culations. 

In the present application where the bare projectile of 
charge Zp collides with a one-electron target of nuclear 
charge Zt, the interaction is given by 

n f z P z T \ 

V{q ' e) = ^9 [~\q + Q(e)/2\ ~ \q-Q(9)/2\) (4) 

where Q(9) can represent any trajectory. Here we use 
a straight-line trajectory to describe the nuclear motion, 
i.e. Q(0) = cosdv + vbsm9b/n, where v is now the 
relative velocity and b is the impact parameter. 
The initial regularized wave function is given by 

t / Q , ^o(v~ 1 q ncsc9 ,-v- 2 ncot9 ) 

MQ, OO) = , „ s , (5) 

X(Q, t>o) 

where ^fo(r,to) is the full wave function at the initial 
time to = — (ficot 9q)/v 2 . At asymptotically large times 
as 9 — > iv, the continuum components of <fr(q, 9) are inde- 
pendent of fl and directly give the ejected electron mo- 
mentum distributions [9. 



A(k) = lim <&(q, 



(6) 



where q — > k/v in the limit and k is the electron wave 
vector. Note that the present formulation does not 
require projection onto dynamic two-center continuum 
states at large times since the bound states shrink to 
vanishingly small regions in q-space. Eq. ^ can be then 
solved using a number of different numerical methods. 

We compute the ejected electron spectra for proton im- 
pact ionization of hydrogen atoms in the impact energy 



range of 5—25 keV where there are long standing contro- 
versies regarding the momentum distribution of ejected 
electrons [14] . Owing to its numerical accuracy and 
efficiency we adopt here the LTDSE technique [51 [71 [5]. 
As the two moving charges become separated by a very 
large distance, the electronic probability density in the 
continuum becomes infinitesimally small, spread over an 
enormous spatial volume. As a result, previous calcula- 
tions have always terminated at relatively small distances 
on the order of 30 to 100 a.u. [7J [51 Q3] before the con- 
tinuum electron density becomes so small that it poses 
intractable numerical difficulties. 

Because of the present singularity-free, scaled repre- 
sentation, the RTDSE does not suffer from this imped- 
iment. Furthermore, due to the steep walls of the har- 
monic oscillator potential in Eq. (3), reflections at the 
lattice boundaries are inherently insignificant. Finally, 
removing the essential singularity also allows us to in- 
tegrate along a path in the complex 9— plane where the 
interaction potential never vanishes for real q, allowing 
us to avoid conventional, ad hoc softcore modifications of 
the potentials near the Coulomb singularities. 

Regarding specific parameters of the present calcula- 
tions, we choose w 1 and the scaled box size < 7, 
i = x, y, z with 256 3 points in order to minimize the level 
of numerical noise in the resulting electron momentum 
distributions. This is judged by performing a test cal- 
culation with projectile charge set to zero. Using the 
quantity 1 — |(<!>(0)|<i>o(0))| 2 as measure of the noise, we 
find it to be less than 3 x 10 -7 . This error is small enough 
so that all of the momentum distribution within a box 
|fc/f| < 2.5 (e.g., see Fig(T| is numerically significant. 
Such small and controllable errors set a standard for ac- 
curacy of full 3-dimensional numerical solutions of the 
TDSE and allow us to see features of electron distribu- 
tions that are beyond the previous state-of-the-art. 

We begin by noting that in this energy range consid- 
ered here, ionization occurs via transfer of energy from 
nuclear to electronic motion in regions where there are 
periodic, unstable classical trajectories. There are two 
types of these trajectories, namely, the superpromotion 
(S) trajectories along the internuclear axis (e.g., as in 
the Fermi Shuttle mechanism [15 ), and the top-of-barrier 
(T) trajectories where the electron moves perpendicular 
to the axis joining the nuclei near the region where the 
net force on the electron vanishes. The S-type motion is 
effective as the nuclei approach while the T-type is effec- 
tive as they recede. Both types of trajectories are related 
to quantal transitions between H^"— like adiabatic molec- 
ular energy curves |16j . 



Figure [TJshows 2-dimensional slices of the ejected electron 



momentum (or wave vector fe) distributions taken with 
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(a) 5 keV, b=0.77, R=52 


(b) 5 keV, b=0.77, R=520 




(c) 5 keV, b=0.77, R=5200 



(d) 5 keV, b=0.77, R=52000 
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FIG. 1: The time-evolution of the spectrum of ejected electrons at a collision energy of 5 keV and an impact parameter of 
&=0.77 a.u. for proton impact ionization of hydrogen. The time integration is terminated at 4 internuclear separations: (a) 
R=52 a.u., (b) i?=520 a.u., (c) _R=5,200 a.u. and (d) i?=52,000 a.u. The electronic amplitude (Eq. Q) has values of ~0.02 
near the edges of the distribution and ~0.64 a.u. at q x =0 and q z ~i0-5. Note that the electron wave vector k is measured in 
units of the relative ion velocity v. 
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FIG. 2: Plot of the ionization amplitude \ A\ of Eq. (JgJ) showing 
target and projectile cusps with the bound states removed. 

the g z -axis parallel to the projectile velocity v, the q x -axis 
parallel to the impact parameter b and y equal to zero, 
computed for collision energy 5 keV with 6=0.77 a.u., and 
with the time integration terminated at four internuclear 
separations: R = 52, 520, 5,200 and 52,000 a.u. One can 
see that the distribution changes significantly between 52 
and 5,200 a.u., predicted in a special limiting case [17] . 
with a small, but noticeable, change between 5,200 and 
52,000 a.u. 

At all post-collisional internuclear distances there 
are fast electrons, i.e. electrons with velocities greater 
than the projectile velocity v, distributed approximately 
isotropically about the center of mass. These electrons 
originate from two S-type transitions near R—0A a.u. 
and R=l a.u. The electrons are distributed in angle ac- 
cording to the angular part of the adiabatic wave 
function at complex values of the internuclear distance 
R [17] . The relevant angular parts have not been re- 
ported in the literature, but they are assumed to be given 



in terms of spherical harmonics with relatively low values 
of angular momentum t. This agrees qualitatively with 
the ab initio distributions reported here. 

The 'elliptical' ring of electrons between the target and 
projectile is another remarkable feature which becomes 
narrower as R increases, and can be interpreted in terms 
of T-type transitions, i.e. electrons are ejected by a pro- 
cess that involves a united-atom rotational coupling tran- 
sition followed by T-promotion to the continuum. It is 
known that united-atom rotational coupling produces a 
state of 7r symmetry with a node along the internuclear 
axis. The shape of this feature at R—52 a.u. and its 
narrowing for increasing R support this interpretation. 

Electron cusps at fe=±u/2 are also seen in Fig. [I] al- 
though one must remember that the bound state elec- 
trons are not separated from the cusps in this figure. It 
is necessary to extrapolate the distribution to small k 
when R is asymptotically large in order to identify these 
features correctly. Figure [2] shows the ionization ampli- 
tude \A\ at i?=52,000 a.u. The electrons in the cusps 
are a small fraction of the total that are ejected, how- 
ever, even this small component can be identified with 
the theoretical approach that we report here for the first 
time with such a time-dependent grid approach. 

The identification of the these features is supported by 
the distributions shown in Fig. [3] for increasing projectile 
velocity. Note that the superpromotion features shrink 
in extent but do not change shape. This is in accord with 
predictions of an exponential exp[— 0A(k/v) 2 v] distribu- 
tion [5] • Alternatively, the T-electron distribution retains 
the same spatial extent but oscillates as identified in Ref. 
[9]. In contrast with that earlier approximate theory, it 
can be seen that some of the electrons actually move from 
the center to the target and projectile nuclei. Note that 
all of the features identified here are superimposed co- 
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(b) 15 keV, b=0.77, R=5200 
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(c) 25 keV, b=0.77, R=5200 
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FIG. 3: Same as in Fig. [T] but for higher collision energies 
with i?=5,200 a.u . The amplitude ranges from ~0.04 near 
the edges of the distributions to ~1.36 a.u. at q x =0 and 
q z =±0.5. 



herently so that interference effects may emerge. Indeed, 
interference structure near the cusps in Fig. [2] is appar- 
ent. A more detailed analysis of the time-dependent wave 
functions shows that the S-type electrons interfere with 
the T-type electrons that have moved from the midpoint 
of the internuclear axis. Such interference has not been 
noted previously since a coherent treatment of S- and 
T-type processes has been beyond the reach of theory. 

The most important conclusion from the present cal- 
culations is that, even for the many-facitcd electron dis- 
tributions seen in Figs. T]|3 in experiment, and in earlier 
calculations [9 , the overall phase factor of the electronic 
wave function in continuum is given by Eq. ^. In {r,t}- 
space this factor for R » 1 is just exp[ir 2 /2t] as derived 
from the free electron propagator. With this phase factor 
removed, a regularized Schrodinger equation is obtained 



that is solved efficiently using the LTDSE method. This 
has allowed us to develop a general theoretical approach 
that allows one to solve the time-dependent Schrodinger 
equation from microscopic to macroscopic distances ac- 
curately. The method is well-suited for treating contin- 
uum electrons in time-dependent fields. Computations 
of ejected electron distributions produced by proton im- 
pact on atomic hydrogen have been used here to illustrate 
the method. The results show measurable features not 
amenable to previous theoretical approaches. The fea- 
tures are in qualitative accord with electron distributions 
measured for proton impact on helium [2]. 
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